Controlling the partial coalescence of a droplet on a vertically vibrated bath 
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A new method is proposed to stop the cascade of partial coalescences of a droplet laid on a liquid 
bath. The strategy consists in vibrating the bath in the vertical direction in order to keep small 
droplets bouncing. Since large droplets are not able to bounce, they partially coalesce until they 
reach a critical size. The system behaves as a low pass filter : droplets smaller than the critical size 
are selected. This size has been investigated as a function of the acceleration and the frequency of 
the bath vibration. Results suggest that the limit size for bouncing is related to the first mode of 
the droplet deformation. 

PACS numbers: 47.55.D-, 47.55. df, 47.85.-g 



I. INTRODUCTION 

Liquid droplets are studied more and more intensively 
in the framework of microfluidic applications [1]. Among 
possible applications, droplets can be used in chemical 
engineering in order to mix tiny amounts of reactive sub- 
stances. To achieve such a goal, one has to invent some 
processes to manipulate a droplet : motion, binary colli- 
sion, fragmentation... To prevent droplets from contam- 
ination, it is also important to avoid any contact with 
a solid element. This paper introduces a method that 
could eventually realize all those operations. The basic 
idea is to combine two physical phenomena : the delayed 
coalescence on a vibrated interface, and the partial coa- 
lescence. 

First, we describe both phenomena separately and 
identify the physical conditions required to achieve these 
effects. 



A. Partial coalescence 

When a droplet is gently laid on a liquid bath at rest, 
it very quickly coalesces with the bath. Because of the 
low air viscosity, the drainage of the air film between the 
droplet and the bath lasts less than a second (a similar 
drainage process can be found in antibubbles Q)- Sud- 
denly, the film becomes thin enough to break and there 
is a contact between the droplet and the bath. A droplet 
made of a low-viscosity liquid can experience a partial co- 
alescence [3|: it does not fully empty, and a new smaller 
droplet is formed above the bath interface (as illustrated 
in the second image line of FiglT]) . This daughter droplet 
can also coalesce partially. The process is repeated un- 
til the droplet becomes small enough to coalesce totally. 
When the droplet is made of a non-viscous liquid sur- 
rounded by air, the daughter droplet has a radius 
that is about half [4] of the mother one (Ri) : 
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[6|. In their experiments, another immiscible liquid was 
surrounding the droplet (instead of air). Recently, many 
studies have focused on understanding this phenomenon 
0, S B H) [S EH ■ In general, gravity and viscosity forces 
in both fluids tend to favor total coalescence Partial 
coalescence is only possible when surface tension is the 
dominant force. A measure of the relative influence of the 
droplet viscosity compared to surface tension is given by 
the Ohnesorge number 
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where v is the viscosity of the liquid, p its density, and 
a the surface tension between the droplet and the air. 
As shown in 0, 0, , the critical Ohnesorge correspond- 
ing to the transition from partial to total coalescence is 
approximately Ohc — 0.026. Small droplets [Oh > Ohc) 
undergo a total coalescence while larger ones [Oh < Ohc) 
coalesce only partially. 
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This phenomenon was first reported in 1930 by Maha- 
jan [5|, and investigated by Charles and Mason in 1960 



FIG. 1: Two droplets arc bouncing on a vibrated liquid/air 
interface. They meet and coalesce. The new droplet is too 
large to take off from the interface and to bounce, so it even- 
tually coalesces partially after 110 ms. The daughter droplet 
is small enough to bounce again for several tens of seconds. 
The six pictures in this sequence are taken every 55 ms. 
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B. Delayed coalescence 

Typically, droplets cannot stay more than a second on 
the liquid /air interface, due to gravity. In 2005, Couder 
et al. |12l . |l3| have found an original way to avoid coales- 
cence. In their experiment, the liquid bath is vertically 
vibrated by an eletromagnetic shaker. The sinusoidal vi- 
bration has a pulsation lo and an amplitude A. When the 
reduced acceleration F = Auj"^ / g is sufficiently large, the 
droplet bounces periodically on the vibrated interface. 
The minimum acceleration Tm for bouncing to occur has 
been estimated by Couder et al. \12\ in the case of viscous 
droplets [1^ as: 
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where Rcc is a critical Reynolds number (constant), pa 
the air density, pi the liquid density, g the gravity accel- 
eration, R the droplet radius (when spherical) and the 
horizontal extension of the film between the droplet and 
the bath. Although rn has no simple analytical form, it 
has been shown ^121] that is monotonically increasing 
with R. This means that for a given reduced acceleration 
r, there is a critical size i?M(r) such as 
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One can reasonably think that this conclusion is still valid 
when droplets are not viscous. This point will be shown 
below. In reality, droplets experience a finite lifetime 
before coalescence, even when they bounce: the film be- 
tween the droplet and the bath breaks and the droplet co- 
alesces sometimes after a thousand bounces. More- 
over, when the vibration amplitude is too high, steady 
waves appear on the bath interface, due to the Faraday 
instabihty [ISj]. Under these conditions, the trajectory of 
the droplet generally becomes chaotic and the lifetime is 
considerably reduced. The more viscous the bath is, the 
higher is the Faraday acceleration threshold F p . Couder 
et al. have noted that when the bath is vibrated just 
below F F , the droplets are sometimes able to move hori- 
zontally on the interface fl3, 16]. In this regime, droplets 
can strongly interact among themselves or with obstacles 

In the present study, the key idea is to combine both 
effects: (i) the partial coalescence, that allows to change 
the radius of the droplet according to Eq. ([T]) and (ii) the 
delayed coalescence that becomes efficient when Ri+i < 
Rm- Conditions on surface tension and viscosity have 
to be satisfied in order to observe both phenomena. Ac- 
cording to the Ohnesorge criterion j24], the kinematic 
viscosity of the droplet liquid has to be lower than 3.7cSt 
in order to ensure the partial coalescence of a millimetric 
droplet (with pi ~ lOOOkg/m^ and cr ~ 20 x lO^^N/m). 
For such a low viscosity, F^^- ^ 1 while F^ > 1 according 
to Eq.®. Therefore, it does not seem possible to stop a 
cascade of partial coalescences when using the same vis- 
cosity for the droplet and for the bath. However, one can 



satisfy Eq.([2]) and Eq.Q by considering a droplet made 
of silicon oil with a low viscosity (0.65cSt) bouncing on 
a bath with a large viscosity (lOOOcSt). 



II. EXPERIMENTAL SET-UP 

We apply vertical vibrations to a 7 mm-thick bath 
of lOOOcSt silicon oil in order to recover the daughter 
droplets of the partial coalescence. The frequency of the 
bath vibration is varied from 40 to lOOHz. The accel- 
eration is tuned from to 2.5g and measured with an 
accelerometer; the precision is estimated below 0.4m/s^, 
g being the gravity. 

During a partial coalescence, the coalescing part of the 
droplet forms a pool at the bath surface. The liquid of 
this pool cannot immediately mix with the surrounding 
highly viscous oil. Consequently, the daughter droplet of 
the partial coalescence has unavoidably to make its firsts 
bounces on the pool formed by its own mother droplet. 
This does not perturbe the bouncing too much. Some- 
times, due to inertial effects, the daughter droplet does 
not succeed in quickly stabilizing its bouncing on the vi- 
brating interface (it was observed even when R < Rm)- 
In this case, another partial coalescence is necessary. 

A fast video recorder (Redlake Motion Pro) is placed 
near the surface. Movies of coalescences have been 
recorded at a rate of up to 1000 frames per second. The 
droplet diameter is measured from the images with an 
absolute error due to the finite pixel size (about 3G/im). 
Since the droplet is constantly deforming, this measure- 
ment is made on the few images where the droplet ap- 
pears the most spherical. 



III. RESULTS 

Figure [T] illustrates a typical scenario observed with 
droplets on a vibrated interface. Two small droplets are 
laid next to each other on the bath. They are small 
enough to take off and bounce. They approach each other 
and coalesce into a large droplet. This resulting droplet 
is now too large to take off, as the thin air film that 
separates it from the bath cannot be regenerated during 
the vibration. The droplet experiences a few oscillations 
before coalescing: it is unstable. The coalescence is par- 
tial and the daughter droplet is small enough to bounce 
again. The system behaves as a low pass filter: partial 
coalescences occur, decreasing the radius of the droplet, 
until the radius of the droplet respects the condition ex- 
pressed in Eq.Q. This experiment is repeated several 
times, with various droplet sizes. 

The radii R of the smallest unstable (open symbols) 
and the largest stable (black symbols) observed droplets 
are plotted as a function of the bath acceleration F in 
FiglH for various frequencies 50, 60, 75 and lOOHz (see 
legend). Between both limit values of R, a mean critical 
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restitution coefficient a of the droplet, such as: 



FIG. 2: Phase diagram (T-R). Critical radius Rm as a func- 
tion of the bath reduced acceleration F = Auj^ /g for various 
frequencies ; 50Hz(») - 60Hz(a) - 75Hz(4) - lOOHz(B). Black 
symbols correspond to the largest observed stable droplets 
while open symbols correspond to the smallest observed un- 
stable droplets. The critical radius 7?m(F,/) is located be- 
tween these two bounds. The dashed vertical line is the min- 
imum acceleration for bouncing. The dashdot horizontal line 
correspond to the asymptotic size Ra{f), the maximum size 
of stable droplets, whatever the acceleration. This size ap- 
proximately correspond to Rr, given by Eq.©. The typical 
error bar is illustrated in the upper left corner of the figure. 



value i?Af (r, /) niay be defined as the maximum value of 
the radius for a stable bouncing droplet. 

The minimum acceleration Tmin = 0.35 needed to sta- 
bilize droplets is lower than 1. It seems impossible to find 
a bouncing droplet on the left of the vertical dashed line. 
The on-off behavior of the stabilization process around 
Tmin bas to bc nuanced. Indeed, the lifetime of a bounc- 
ing droplet should be considered close to the boundaries. 
Some combinations of the parameters probably stabilize 
droplets for a longer time than others, or allow to stop the 
cascade of partial coalescences at an intermediate value 
between Rm = and Rm = 0.4mm. 
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This value is higher than the restitution coefficient of 0.22 
found in previous works [20l [2]|. However, some aspects 
have to be taken into account : (1) In those experiments, 
small (less deformable) droplets bounce on a low-viscous 
(deformable) bath. The restitution mechanism is conse- 
quently not the same as in our case. (2) Moreover, for 
elastic beads, it is known that the restitution coefficient 
tends to 1 when the impact velocity tends to [l^]. In 
our case, impact velocities near the critical acceleration 
are about v = Tg/uj < 2cm/s while in (2l| . impact ve- 
locities are at least ten times higher. This could also 
explain the difference between restitution coefficients in 
both configurations. 

Second, experimental points of Fig[2]define two regions 
in the T — R diagram. The boundaries that depend on 
the frequency are shown as continuous curves based on 
a qualitative fit. For a given frequency, droplets located 
above the curve are unstable while droplets located below 
are stable. More precisely, when a droplet is created with 
a radius larger than the critical droplet radius i?A/(F, /), 
the droplet partially coalesces until it reaches the sta- 
ble region (below the curves). The stability threshold 
RiviiX: I) slightly increases with the forcing acceleration, 
as in Eq.Q. The curves seem to saturate at high accel- 
erations towards an asymptotical value Ra{f) that de- 
creases with an increasing frequency. 

Since the bath is much more viscous than the droplet, 
the droplet deformation is the main elastic mechanism 
able to ensure the bouncing. A characteristic size related 
to the deformation is given by the wavelength of the first 
normal mode. According to for a sitting droplet, 
this wavelength A satisfies to: 
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IV. DISCUSSION 

In our experiments, low viscous droplets are bouncing 
on a high viscous bath. It means that the deformation 
of the droplet is larger than the deformation of the bath 
: the droplet is more able to store and restore interfacial 
energy in order to ensure the bouncing. Conditions are 
comparable to an elastic ball bouncing on a rigid plate. 

First, according to Eq.jH]), Tmin = 1 for viscous 
droplets. However, since non- viscous droplets can restore 
surface energy (the dissipation being limited), the bounc- 
ing is possible for reduced accelerations lower than 1: A 
restitution coefficient has to be taken into account. The 
analogy between the non-viscous droplet and the elastic 
ball (defined and analyzed in [l^) allows us to define this 



The radius Rr corresponding to A is represented in FiglSl 
It fits well to the asymptotical value Ra for high acceler- 
ations. This fact suggests that the first mode of droplets 
deformation is related to the maximum bouncing droplet 
size for a given frequency. 

Finally, the critical radius Rm has been measured as 
a function of the frequency, the acceleration being fixed 
F = 1 (Figl3]). Open circles correspond to the small- 
est unstable droplets, while black circles are related to 
the largest stable ones. The continuous line represents 
Eq.®. The curve does not fit the experimental re- 
sults because the asymptotic regime is not reached when 
F = 1. Moreover, sitting and bouncing droplet geome- 
tries are not exactly equivalent. However, we found 
that the decreasing of Rm with increasing frequencies 
(at F = 1) is coherent with the decreasing of Rr- 
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FIG. 3: Critical radius Rm as a function of the frequency 
when r = 1. Black circles correspond to the largest observed 
stable droplets while open circles correspond to the smallest 
observed unstable droplets. The critical radius Rm is located 
between these two bounds. The continuous line corresponds 
to Eq.®. 



V. CONCLUSIONS 



In summary, we have developed a method to stop a 
cascade of partial coalescences and to maintain the resid- 
ual droplet alive on the interface. Our approach is based 
on the experiment of Couder that consists in the 
bouncing of a droplet on a bath subjected to a vertical 
vibration. We have shown that it is essential to use a 
viscous liquid for the bath and a low-viscous one for the 
droplet. The threshold radius below which the droplet 
can bounce has been found to depend on the accelera- 
tion and on the frequency of the bath. By tuning these 
forcing parameters, it is possible to choose the maximum 
allowed size for bouncing droplets in a broad range. A 
concordance is observed between the radius of the largest 
stable droplet (for a given frequency) and the radius re- 
lated to the first normal mode of deformation. 
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